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PREFACE 


This  Memorandum  Is  a  part  of  continuing  RAND  studies  on  the  be¬ 
havior  of  bodies  re-entering  the  atmosphere  at  high  velocities.  The 
results  should  be  of  Interest  to  those  engaged  In  predicting,  under¬ 
standing,  and  evaluating  the  characteristics  of  hypersonic  trails  In 
the  atmosphere. 


SUMMARY 


An  attempt  Is  made  to  develop  a  theory  for  the  calculation  of 
the  growth  of  the  turbulent  wake  formed  behind  a  blunt  body  moving 
at  hypersonic  speeds.  The  classical,  semlemplrlcal  subsonic -wake 
theory  based  on  the  notion  of  eddy  dlffuslvlty  Is  used  throughout, 
duly  modified  to  account  for  compressibility  effects.  The  gas  Is 
assumed  to  be  In  thermodynamic  equilibrium.  A  number  of  suitable 
approximations  are  made,  resulting  In  a  closed-form  solution. 

Several  plots  showing  the  growth  of  the  width  of  the  turbulent 
core  downstream  are  presented  for  different  Mach  numbers  and  differ¬ 
ent  values  of  the  enthalpy  prevailing  at  the  center  of  the  "neck" 
of  the  wake.  It  Is  found  that  compressibility  effects  for  blunt 
bodies  are  limited  In  general  up  to  about  1000  diameters  downstream. 
Beyond  this  distance  the  wake  grows  with  the  1/3  power  rule,  like  a 
subsonic  Incompressible  wake,  with  no  dependence  on  the  Mach  number. 

The  fundamental  equations  are  also  written  for  sharp  bodies  for 
which  the  boundary- layer  effects  are  Important.  It  Is  found  that  the 
enthalpy  decay  along  the  axis  downstream  is  faster  for  thin  bodies 
than  It  is  for  blunted  ones. 
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SYMBOLS 

A  >  excess  amplitude  of  turbulent -enthalpy  profile,  defined  In 
Eq.  (19) 

a  ■  mnaerlcal  constant  used  In  Eqs.  (2)  and  (10) 

B  ■  parabolic  "width"  of  turbulent-enthalpy  profile,  defined 
In  Eq.  (20) 

C|^  *■  drag  coefficient 
d  body  diameter 

f  ■  velocity  ratio  defined  In  Eq.  (5) 

H  ■  nondimens lonal  enthalpy,  defined  In  Eq.  (7) 
h  enthalpy 

hf  ■  hj/h^ 

ft 

M  free -stream  Mach  number 

p  ■  pressure 

R  <■  Howarth  variable  In  the  radial  direction,  defined  in  Eq.  (4) 

R^  ■  nondimens lonal  co-ordinate  describing  bow  shock  wave 

■  turbulent  variable  In  the  radial  direction,  defined  In  Eq.  (14) 

r  ■>  radial  co-ordinate 

r_  ■  nose  radius 

o 

r  -  r/r„ 

o 

■  free-stream  velocity 

u  ■■  velocity  In  the  x  direction 

X  ■  co-ordinate  along  the  direction  of  the  free  stream 

■  distance  downstreasi  at  which  the  pressure  is  aad>lent  for 

every  streamline,  given  by  Eq.  (1) 

X  ■  x/d 

Y  “  ratio  of  specific  heats,  c^/c^ 
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"  numerical  constant  used  In  Eq.  (2) 

6  *■  compressible  eddy  dlffuslvlty  defined  In  Eq.  (12) 
e  ■  eddy  dlffuslvlty  defined  In  Eq.  (15) 

.<  ■  numerical  coefficient  used  for  the  definition  of  the  eddy 
dlffuslvlty 

X  ■■  fraction  of  the  stagnation  enthalpy  calculated  at  the  "fieck" 
of  the  turbulent  core;  see  Eq.  (35) 

p  ■  mass  density 
SUBSCRIPTS 

f  ■  the  turbulent  front 
1  ■  an  Initial  condition 
“  ■  the  free  stream 
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I.  lOTBDDUCTION 

It  Is  a  matter  of  conmon  observation  that  the  motion  of  bodies 
In  fluids  produces  a  trail  which  very  often  persists  over  distances 
many  times  the  characteristic  length  of  the  body.  Theoretical  and 
experimental  studies  consist  primarily  of  descriptions  of  the  distri¬ 
bution  of  momentum  defect  downstream.  At  low  speeds  this  defect  Is 
associated  with  the  viscous  layer  only.  At  supersonic  speeds,  how¬ 
ever,  the  major  part  of  the  momentum  defect  for  a  blunt  object  Is  due 
to  the  Irreversibility  created  by  the  shock  wave.  The  maximum  Irre¬ 
versibility  (and  hence  the  maximum  defect)  for  an  axlsymmetrlc  body 
appears  on  the  axis  parallel  to  the  flow  where  the  shock  Is  normal, 
with  smaller  Irreversibilities  and  hence  smaller  momentum  defects, 
corresponding  to  streamlines  located  farther  away  from  the  axis. 

The  structure  of  the  flow  for  a  spherical  object  moving  at  hyper¬ 
sonic  speeds  can  be  seen  In  Fig.  1.  The  fluid  Is  compressed  through 
the  bow  shock  tmve  and  expanded  Isentroplcally  until  It  Is  recompressed 
again  by  the  so-called  trailing  shock.  Then  It  expands  again  until 
all  streamlines  reach  the  ambient  pressure  at  a  distance  denoted  by 
X  ■  x^.  This  distance  can  be  calculated  as  Indicated  In  Ref.  1  by 
using  the  blast-wave  theory,  although  this  theory  allows  only  for  the 
bow  shock  wave  and  neglects  the  details  of  the  flow  from  the  bow  shock 
to  the  point  x  ■  x^.  The  latter  has  been  shown  In  Refs.  1  and  2*  to 
create  a  small  amount  of  Irreversibility  when  compared  with  the  bow 
shock  wave,  so  that  If  one  Is  Interested  In  calculating  the  flow  at 
a  distance  far  enough  dotmstream.  Its  Influence  can  be  neglected  In 
a  first  approximation.  In  nondlmenslonal  terms,  the  distance  x^  Is 
given  In  Ref.  1  as 


d 
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(1) 


In  the  above,  d  Is  the  dlaiseter  of  the  sphere  and  M  Is  the  free-stream 


Also  private  comaunlcatlon  from  L.  Hromas,  Space  Technology 
Laboratories. 


ake  behind  a  sphere 
eeds 
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Mach  ntjmber.  Feldman,  in  Ref.  2,  was  the  first  to  analyze  the  flow 
of  a  hypersonic  trail  around  a  blunt  body.  By  standard  numerical 
techniques  he  calculated  the  flow  around  a  hemisphere  followed  by  a 
cylinder  by  assuming  inviscid  flow  of  a  real  gas  in  thermodynamic 
equilibrium.  At  the  cross  section  x  ■  he  assumed  that  the  enthalpy 
profile  is  given  by  a  Gaussian  distribution,  the  width  of  i^ich  he  ob> 
tained  by  a  best  fit  of  his  numerical  results.  From  there  on  he  as¬ 
sumed  that  further  cooling  occurred  by  the  mechanism  of  heat  conduc¬ 
tion  only.  He  was  thus  able  to  determine  all  thermodynamic  parameters. 
Including  electron  concentration,  up  to  the  point  vdiere  the  ambient 
state  is  reached. 

In  Ref.  1  it  was  shovm  that  the  inviscid  enthalpy  profile  at  the 
cross  section  x  ■  x^  can  be  approximated  analytically  in  a  closed  form 
by  assuming  a  parabolic  bow  shock  wave  as  given  by  the  cylindrical- 
blast  -wave  theory.  The  result  is 

(2) 


In  the  above,  R^  is  the  nondimens ional  radius  (with  respect  to  the 
radius  of  the  nose  r^)  determining  a  stream  tube  entering  the  bow 
shock  wave  (see  Fig.  1)  and  is  related  to  the  nondimens ional  radius 
r  *  r/r^  of  the  same  stream  tube  at  the  cross  section  x  ■  x  by  the 
equation  of  mass  continuity 

P_U  R  dR  ■  P  u  rdr  (3) 

•®S8  XX 

o  o 

Feldman  Indicated  from  his  numerical  calculations  that  the  ratio 

u^  ^  0.8,  so  in  this  sense  the  radius  R  is  almost  equivalent  to 
o 

the  Howarth  variable  R,  given  as 


(4) 
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where 

f  ■  u  AJ,  (5) 

’‘o 

In  Appendix  A  It  ie  shotm  that  £  Is  given  apprexlaately  In  a 
clesed  form  by  the  expression 


For  Mach  numbers  between  IS  and  3S>  f  varies  between  0.77  and  0.83, 
and  It  Is  thus  seen  that  the  Mach<nvumber  dependence  Is  Indeed  very 
weak.  From  the  condition  that  the  stagnation  enthalpy  remains  con* 
stant  everyidiere  In  the  Invlscld  flow  field,  ve  obtain 


For  the  quantity  £,  Eq.  (6)  could  be  substituted  above,  thus  obtain* 
Ing  Eq.  (A*2)  again;  however.  If  we  enter  the  value  £  *■  0.8,  the  re* 
suit  Is  a  good  approximation  for  Mach  numbers  of  10  and  above  to  the 
exact  values  that  can  be  obtained  by  using  real-gas  effects.  Also, 
by  neglecting  unity  compared  to  the  term  containing  the  Mach  number 
In  Eq.  (7),  we  obtain 

»  “  (w 

The  parameter  a  In  Eq.  (2)  and  the  drag  coefficient  Cj^  determine  the 
shape  of  the  shock  wave  from  blast-wave  theory  to  be 

».  -  m 

o 

This  coefficient  accounts  only  for  the  pressure  drag. 
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The  parameter  Y.  Is  an  effective  c  /c  ratio,  which  as  shown  In  Ref.  1 

L  P  V 

has  the  approximate  value  of  1.2.  The  value  of  a  Is  approximately 
equal  to  2.6. 

In  Ref.  1  It  was  also  shown  that  a  good  approximation  of  Eq.  (2) 
Is 


h  (R) 

Jia _ 


In  the  above,  the  values  of  f  0.8,  a  ■  2.6,  and  Y^^  *  1.2  have  been 
used.  This  relation  Is  In  good  agreement  with  numerical  results  ob¬ 
tained  through  the  method  of  characteristics  In  Refs.  2  and  3.  A 
comparison  Is  shown  In  Fig.  2,  where  the  exact  value  for  H  was  used 
In  Eq.  (10)  rather  than  the  approximation  of  Eq.  (8). 

Further  Inspection  of  Fig.  1  shows  that  the  trailing  shock  Is 
formed  around  a  "neck”  which  presumably  contains  the  mass  flowing 
through  the  boundary  layer  over  the  sphere.  From  the  point  of  view 
of  total  momentum  defect  (and  hence  drag)  contained  In  the  neck,  sim¬ 
ple  calculations  by  Lees  and  Hromas  show  that  for  a  typical  hypersonic 

(3) 

flight  this  amount  Is  very  small  compared  with  the  total  drag. 

They  point  out  that  this  drag  represents  the  sum  of  the  skln-frlctlon 
drag  on  the  body  and  the  momentum  defect  associated  with  the  pressure 
rise  across  the  trailing  shock  at  the  neck.  Their  estimate  of  these 
two  contributions  at  the  neck  Is  given  for  the  cylindrical  case  as 


H 


1  + 


•Vb 


(10) 


^^D^neck 


(11) 


This  result  should  be  multiplied  by  for  the  axlsymsetric  case. 
The  Reynolds  number  Is  based  on  the  dlaswter  and  free-stream  condi¬ 
tions  . 

Turbulence  Is  always  associated  with  regions  of  high  vortlclty, 
which  In  the  present  case  Is  the  free  shear  layer  formed  after 
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Fig.  2  —  Inviscid  enthalpy  profiles  from  Ref  3  compared  with  present 

work  for  M  =  8.5 


7 


separation  occurs.  Although  conclusive  experimental  study  of  the  so- 
called  "near  wake"  region  Is  still  lacking,*  we  shall  assume  that 
turbulence  Is  present  In  the  region  near  the  neck,  engulfing  the  fluid 
mass  downstream  that  has  gone  over  the  blunt  body  through  the  Invlscld 
part  of  the  wake.  The  "engulfing"  mechanism  for  the  first  few  diam¬ 
eters  downstream  of  the  neck  will  be  dictated  solely  by  expansion, 
since  Just  behind  the  trailing  shock  the  temperature  and  pressure  are 
very  high.  As  the  pressure  drops  to  nearly  ambient,  turbulent  diffu¬ 
sion  will  become  predominant— swallowing  up  more  and  more  momentum  and 
finally  spreading  to  a  point  where  the  whole  momentum  deficiency  Is 
associated  with  turbulence. 

A  calculation  based  on  the  above  model  needs  the  following  g. 
priori  knowledge;  (a)  the  position  and  thickness  (In  terms  of  mass, 
momentum,  and  energy)  of  the  neck,  (b)  pressure  profiles  downstream 
of  the  neck,  and  (c)  a  model  for  the  estimation  of  the  rate  of  energy 
diffusion  due  to  turbulence.  The  first  two  Items  presume  the  solution 
of  the  near-wake  problem  with  an  answer  describing  the  Interaction  of 
the  trailing  shock  with  the  free  shear  layer,  a  rather  difficult  and 
still  unsolved  problem.  As  a  rough  approximation  one  may  obtain  some 
Information  for  the  first  two  Items  by  solving  the  Invlscld  problem 
over  a  sphere  followed  by  a  cone,  which  takes  the  place  of  the  shear 
layer,  and  then  a  cylinder  fitted  at  the  neck.  Shadowgraph  pictures 
can  be  used  for  guidance  In  prescribing  the  above  geometry.  Such  cal¬ 
culations  are  reported  In  Ref.  3.  For  Item  (c)  the  only  choice  of¬ 
fered  Is  the  extension  of  the  subsonic  eddy-dlf fuslvlty  theories  to 
the  compressible  case.  Such  an  extension  Is  provided  In  Refs.  3  and 
5.  Both  references  present  an  equivalent  eddy  dlffuslvlty  $  given  In 
terms  of  the  Incompressible  eddy  dlffuslvlty  e  as 


2 


Some  preliminary  results  are  reported  In  Ref.  4. 
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In  the  above,  P(0)  indicates  the  nass  density  at  the  axis,  and  the 
density  at  the  "front"  of  the  turbulent  core  (see  Fig.  1);  j  takes 
the  value  0  for  the  t%K)-dimsnsianal  case  and  1  for  the  axisymatric. 

At  this  point  let  it  be  stated  that  there  is  no  rigorous  theore¬ 
tical  understanding  even  of  the  rudiaents  of  shear  turbulence;  uhat 
has  been  accomplished  over  so  many  years  of  study  is  the  replacement 
of  our  ignorance  by  one  or  more  "factors  of  adjustmnt"  which  can  be 
conveniently  selected  to  fit  the  experiments.  For  instance,  for  tur¬ 
bulent  flow  in  a  pipe  we  have  the  numerical  coefficient  of  proportion¬ 
ality  yielding  the  linear  law  of  variation  of  the  mixing  length, 
tihereas  for  wakes  we  select  the  numerical  value  of  a  constant  of  pro¬ 
portionality  yielding  e,  a  number  akin  to  a  Reynolds  number.  Mow  e 
(as  will  be  seen  later)  is  also  proportional  to  the  drag  coefficient, 
and  in  the  problem  at  hand,  the  drag  associated  with  the  turbulent 
core  is  variable  in  the  direction  of  the  stream.  It  is  therefore 
suggested  in  Ref.  3  that  experimental  evidence  should  also  be  used 
in  deciding  between  the  two  extreme  cases  of  a  coefficient  propor¬ 
tional  to  a  frozen  drag  (say  the  oi\e  at  the  neck)  or  local  drag. 

From  a  logical  point  of  view  one  necessarily  concludes  that  the 
artificial  mechanisms  of  eddy  diffusivlty  and  mixing  length  explain 
nothing  of  the  true  nature  of  turbulence.  In  essence  they  are  used 
to  provide  definitions  and  numerical  values  for  quantities  whose  in¬ 
trinsic  existence  can  be  justified  only  by  the  fact  that  they  some¬ 
times  provide  adequate  answers  in  much  wider  regions  of  observation 
without  additional  forcing,  although  they  are  adjusted  to  satisfy  the 
experimental  evidence  in  a  restricted  area.  From  this  point  of  view, 
the  merit  of  a  scmiempirlcal  theory  is  judged  by  the  degree  of  uni¬ 
versality  of  the  answer  for  a  minimum  number  of  factors  of  adjustment. 

Bearing  this  in  mind,  we  will  attack  the  problem  at  hand  by  as¬ 
suming  that  the  inviscid-enthalpy  profile  prevailing  in  the  "laminar 
region"  is  given  by  Eq.  (10).  Inside  the  turbulent  core  we  assume  a 
parabolic  profile  extending  up  to  the  distance  r^  characterizing  the 
turbulent  front. 


In  contrast  to  the  subsonic  case  in  which  the  turbulent  core 
encompasses  all  the  drag. 
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In  the  limit  r  -*  0,  iihere  the  energy  equation  will  be  satisfied 
in  detail  along  the  axis,  a  parabolic  profile  is  the  only  one  satis¬ 
fying  this  equation.  In  the  radial  direction,  energy  will  be  con¬ 
served  in  bulk  between  r  ■  0  and  r  ■  in  the  fashion  of  the  well- 
known  boundary-layer  integral  methods.  We  will  compute  eddy  dlffu- 
slvity  through  Eq.  (12),  using  for  e  its  incompressible  definition, 
and  bearing  in  mind  that  this  definition  is  strictly  valid  only  in 
the  neighborhood  of  the  axis  of  symmetry. 

In  every  physical  detail  the  above  is  the  physical  model  adopted 
(3) 

by  Lees  and  Hromas;  ho%iever,  the  mathematical  model  and  approxi¬ 
mations  of  the  present  %K>rk  are  different,  yielding  a  closed-form 
solution  rather  than  a  solution  necessitating  extensive  numerical 
calculations. 
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III  ABP  IBL  gAg&  er.  VUmS  BOpiBg 

The  equetlen  of  eeneervetlen  of  energy  for  an  axlaymetrle  wake 
is 

Here  again  e  Is  the  Bousslnesq  turbulent  eddy  dlffuslvlty.  For  the 
velocity  Ut  we  set 


OD 


where  f  is  between  0.8  and  1.0. 

In  Eq.  (13)  we  have  neglected  the  contribution  of  the  nomal  ve¬ 
locity  component  in  the  convective  heat  transfer*  any  cooling  due  to 
pure  expansion*  and  also  the  contribution  due  to  molecular  heat  con¬ 
duction.  We  now  define  the  turbulent  Howarth  variable  from  the 
transformation 


(14) 


In  the  above,  p.  is  the  mass  density  prevailing  at  the  turbulent  front. 

*  (6) 

From  the  semiempirical  wake  theory*  we  set^  ' 


€•  K 


(15) 


From  the  approximate  conservation  of  the  stagnation  enthalpy*  we  also 
have 


h  -  h^ 


-  u^  (U,  -  u)(U,  +  u) 
2  "  2 


-  »> 


(16) 
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Writing  the  ab«ve  at  r  ■  0  and  r  ■  we  have 


h(0)  -  hj 


Substitution  In  Eq.  (15)  yields  (for  J  *  1) 


«  No)  -I’f1 


(17) 


(18) 


Now  let  the  turbulent  profile,  assumed  to  be  parabolic,  be  expressed 
as  an  excess  of  the  enthalpy  at  the  turbulent  front  by  the  equation 


where 


In  the  neighborhood  of  the  axis  we  can  set 

2 

-  Pr 


(19) 


(20) 


(21) 


We  new  use  Eqs.  (19)  and  (21)  In  Eq.  (18)  to  calculate  6  at  the  axis, 
since  the  equation  of  energy  conservation  will  be  satisfied  In  detail 
only  there.  The  result  Is 


U  d 

OD 


A(x)rj 


(H^)ui[pfJ 


(22) 
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The  above  was  first  derived  In  Ref.  3.  At  the  axis,  Eq.  (13)  yields, 
after  introduction  of  the  variable  R^,  the  diffuslvlty  ^ ,  and  the 
approximation  of  Eq.  (21),  the  following: 


Sx 


fU^d 


w  ** 


e 


(23) 


The  factor  4  In  Eq.  (23)  enters  because  x  la  nendlnenslonallaed  with 
respect  to  the  nose  diameter  d,  whereas  R  Is  nandlmenslonallzed  with 
respect  to  r^. 

The  quantity  In  brackets  above  can  be  calculated  easily  from 
Eqs.  (19)  and  (20).  The  result  Is 

^  .  -aA(x)  g 

In  the  neighborhood  of  r  •  0  we  have 


(24) 


(25) 


and  hence 


fVd  ,2  ® 


(26) 


The  left-hand  side  of  this  equation  is  now  calculated  at  r  ■  0  with 
the  help  of  Eq.  (19) 


As  shown  In  Ref.  3  from  a  conservatlon-of^nass  argument,  we 

have  Rj  -  (P-/P,)l4  -  <Y  -  DH^C  /8(1  +  H).  If  we  use  Eq.  (8), 
f  fi 

the  quantity  In  the  right-hand  side  Is  approximately  eq\ial  to 
(3/4)Cjjj^,  a  very  small  quantity  compared  with  the  other  terms  except 

very  close  to  the  neck. 
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aix 


(27) 


Appendix  B  shows  that  It  Is  legltlnate  to  neglect  the  reduction  of 
the  Invlscld  enthalpy  at  the  front  due  to  the  siechanlso  of  molecular 
conductloui  so  that  from  Eq.  (10)  we  obtain 


% 

dx 


aA(R ) 


(29) 


Before  we  can  solve  this  equation  for  ve  need  to  evaluate 
the  function  A(R^).  This  Is  done  In  Appendix  C,  where  a  calculation 
la  Blade  of  the  amount  of  energy  (or  equivalent  drag  In  terms  of  the 
local  drag  coefficient  )  contained  within  the  turbulent  core  at 
a  given  station  x.  The  tital  aaiount  of  energy  (or  total  amount  of 
drag)  between  r  ■  0  and  r  "  remains,  of  course,  constant  along  x, 
so  that  the  difference  In  the  drag  of  the  turbulent  core  at  t%K>  dif¬ 
ferent  axial  stations  Is  exactly  equal  to  the  difference  of  the  energy 

(or  drag)  contained  by  the  Invlscld-energy  profiles  at  the  same  two 
* 

stations.  After  some  approximations  In  Appendix  C,  we  obtain 
Eq.  (C-7),  which  for  convenience  Is  combined  here  with  Eq.  (C-4). 


A(R^0 


C  (Y  -  1)H^ 

Jll _ 


Cv  -  1)H^ 
r2  R^ 

fl 


(29) 


We  now  use  In  Eq.  (24)  the  expressions  for  A(R^)  and  6  from  Eqs.  (29) 
and  (22).  After  soise  algebra  we  find 


Uniform  pressure  la  assumed  In  this  argument. 


This  equation  gives  the  growth  of  the  turbulent  core  In  the  Howarth 


In  this  step  we  add  one  to  the  right-hand  side  of  Eq.  (10)  in 
order  to  make  it  an  exact  eqiiatlon  as  R  ^  since  the  density  ratio 
nust  be  accurate  for  both  small  and  large  R. 


IS 


plane.  In  Appendix  by  using  the  Ideal  gas  law^  we  calculate 
as  a  function  of  the  physical  co-ordinate  after  using  Eq.  (10). 
The  result (  repeated  here.  Is 


Now  the  function  r,(x)  could  be  calculated  If  we  knew  the  quantities 

*  2 
p(0)/p^,  the  constant  of  Integration, <  ,  and  /C^  . 

As  has  been  mentioned  before,  the  pressure  ratio  p(0)/p^  can 
only  be  given  by  solving  exactly  the  Interaction  problem  of  the  trail¬ 
ing  shock  with  the  wake  core.  Since  such  a  solution  Is  not  available, 
the  approximations  of  Ref.  3  described  In  the  Introduction  of  this 
Nemorandum  are  appropriate.  The  constant  of  Integration  appearing  In 
Eq.  (32)  can  be  set  equal  to  a  value  giving  the  position  of  the  neck 
and  Its  thickness  as  observed  In  experiments;  this  was  also  the  pro¬ 
cedure  adopted  In  Ref.  3.  It  Is  obvious  that  the  value  of  this  con¬ 
stant  has  an  Influence  only  a  very  few  diameters  downstream  from  the 
neck. 

The  value  of  k  Is  again  offered  as  an  additional  arbitrary  param¬ 
eter,  but  the  remarkable  thing  Is  (as  mentioned  In  Ref.  3)  that  ex¬ 
perimental  evidence  suggests  the  same  value  used  In  subsonic  wakes. 
This  value  Is  approximately  equal  to  0.035. 

2 

Finally,  an  estimate  must  be  made  of  the  quantity  R,.  /C_  . 

ri  Ugj 

From  Eq.  (29),  recalling  the  definition  for  A(x)  from  Eq.  (19),  we 
find 


For  the  ratio  h../h^  we  my  use  the  value  h  /h^,  given  by  Eq.  (10), 

*o 

since  Initially  we  are  very  close  to  the  axis. 
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An  exact  cmiputatlen  of  the  ratio  h^(0)/h^  Is  not  possible »  since 
the  near-wake  preblem  remains  unsolved.  However,  the  physical  de¬ 
scription  of  Bef.  7  could  be  used  as  a  guide.  If  we  assume  that  the 
"dividing  streamline"  brings  the  fluid  to  rest  Just  at  the  center  of 
the  neck,  the  enthalpy  for  an  Insulated  body  will  be  approximately 
equal  to  the  stagnation  enthalpy  (for  Frandtl  number  of  order  one). 

If  heat  transfer  Is  allowed  at  the  body,  then  this  enthalpy  will  be 
smaller. 

it 

Now  let  X  be  the  fraction  of  the  stagnation  enthalpy  to  which 

h. (0)/h^  Is  equal.  For  the  reason  explained  above  we  shall  treat  X 

2 

as  an  additional  free  parameter.  The  ratio  R,.  /C.  depends  on  X 
only  and  can  be  calculated  from  Eq.  (34)  after  use  Is  made  of  Eq.  (8): 


""fl 


(35) 


2 

For  the  values  of  R  /C_  ■  3,  4,  6,  8,  10,  and  12,  the  corre- 

tl 

spending  X's  ■  1.0,  0.83,  0.67,  0.58,  0.53,  and  0.50. 

Ue  proceed  now  to  examine  the  asymptotic  character  of  Eq.  (32) 
for  small  and  large  distances  downstream.  It  Is  easy  to  establish 
that  for  small  values  of  R^,  Eq.  (32)  gives 


R 


3 

f 


dx  +  (const) 


(36) 


If  we  approximate  Eq.  (33)  for  small  R^'s  and  use  Eqs.  (8)  and  (35), 
the  above  becoises 


(const)  (37) 


The  stagnation  enthalpy  Is  equal  to  h^^/h^  (y  -  l)H^/2  3H. 
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At  high  values  of  irfien  the  influence 
ceases  to  exist,  Eq.  (32)  becemes 


of  the  parameter 


2 


r 


3 

f 


UKC^ 


(38) 


2 

It  is  thus  concluded  that  the  value  of  R,.  /C.  and  H  (and  hence 

ri 

the  Mach  number)  is  only  significant  during  the  initial  period  of  the 

process.  Very  far  downstream,  only  the  drag  (essentially  equal  to 

the  total  drag,  since  C-  is  so  small)  affects  the  growth  of  the 

“fi 

wake,  exactly  as  in  the  subsonic  case.  For  all  the  cases  studied, 

this  distance  is  of  the  order  of  1000  diameters  downstream.  In  a 

log-log  presentation,  the  straight  line  of  Eq.  (37)  can  assume  a 

higher  or  lower  position  relative  to  Eq.  (38).  The  position  will  be 

higher  at  the  higher  Mach  numbers  and  higher  values  of  \  (no  heat- 

transfer  effects),  as  one  would  expect  on  physical  grounds. 

The  function  r.(x,  R,.^/C-  )  for  M  •  8.5  is  plotted  in  Fig.  3 

fl 

for  the  several  values  of  X,  along  with  the  experimental  data  of 

Ref.  8.  It  is  seen  that  the  Influence  of  the  factor  X  is  limited 

at  high  X,  and  certainly  an  exact  value  for  X  or  R,.  /C-  cannot  be 

£1 

decided  from  the  experiments,  since  the  scatter  is  high  at  the  smaller 

x's  where  the  Influence  of  X  is  stronger.  For  Fig.  3,  the  constant 

of  integration  of  Eq.  (32)  was  set  equal  to  zero.  Given  a  set  of 

values  of  C-  and  X,  the  quantity  R,.  is  fixed,  and  hence  so  is  r,. . 
fl 

For  different  values  of  C.  but  fixed  X's  we  will  have  different 

r^^'s,  which,  however,  should  physically  correspond  to  the  same  x^. 

The  value  of  x^  should  be  selected  from  experimental  observation,  and 
the  constant  of  integration  accordingly  fixed  so  that  the  curve  goes 
through  the  point  r^^,  x^.  This  procedure  amounts  to  an  appropriate 
sliding  of  the  co-ordinate  x  by  a  constant  (positive  or  negative) 
value.  In  Fig.  3  the  curve  obtained  in  Ref.  3  allowing  for  pressure 
variation  is  shown,  along  with  another  one  in  which  the  pressure  was 


001 


L8 


Fig.  3  —  Comparison  of  the  growth  of  the  turbulent  core  of  o  blunt  body  with  Cq-I.O  and 
several  values  according  to  present  theory  (constant  pressure)  with  theory  of  Ref.  3 

and  the  experimental  doto  of  Hef.  8  for  M-8.5 


19 


kept  constant  everywhere.  It  should  be  noted  that  the  curves  of 
Fig.  3  obtained  in  Ref.  3  correspond  to  the  value  R  ^/C_  “•  10; 

II 

It  Is  seen  that  the  agreement  with  Eq.  (32)  is  very  good. 

In  Fig.  4  the  function  r,(x)  Is  shown  plotted  for  different 
2  * 

values  for  R,.  /C_  and  M  ■  22.  Here  again  the  constant  of  Inte- 
fl 

gratlon  In  Eq.  (32)  Is  set  equal  to  zero.  In  the  same  figure  the 

(3) 

curves  calculated  for  constant  and  variable  pressure  are  shown. 

2 

The  quantity  R,.  /C_  was  taken  in  Ref.  3  to  be  equal  to  about  11. 
fi 

It  Is  apparent  that  even  when  adjustment  Is  made  for  the  curve  cor- 

responding  to  R,.  /C_  ■  10  or  12  to  pass  through  the  point  x  ■  7.5 

11 

and  r^  ■  1.8,  the  agreement  with  Ref.  3  Is  not  as  good  In  the  present 
case  as  It  was  for  the  case  M  ■  8.5.  It  should  be  noted,  hotfever, 
that  the  general  trend  of  the  curve  r^(x)  is  the  same  as  the  one  of 
Ref.  3.  The  numerical  differences  (certainly  small  when  considering 
the  scatter  of  the  usual  experimental  data  obtained  In  the  ballistic 
range  and  the  crude  character  of  the  semiemplrlcal  turbulent  theory 
used)  are  due  to  the  fact  that  the  function  r£(x)  depends  on  the 
slopes  of  both  functions  h^  and  A,  as  can  be  seen  from  Eq.  (27). 

Since  both  h^  and  A  (Eqs.  (10)  and  (29))  provide  inflection  points, 
differences  are  to  be  expected  when  coiaparlng  results  with  other  ana¬ 
lytical  representations  of  these  curves,  which  are  in  good  agreement 
in  the  ordinates  but  have  significant  differences  in  the  slopes. 

In  Fig.  5  a  plot  similar  to  Figs.  3  and  4  is  shown  for  M  •  16. 
Observation  of  Figs.  3,  4,  and  5  shows  that  for  small  and  very 
high  X  the  turbulent  core  (pressure  variation  apart)  grows  with  the 
1/3  power  of  x  according  to  Eqs.  (37)  and  (38).  In  Fig.  6  the  value 
of  the  quantity  R,.  /C-  is  kept  constant  (equal  to  10),  and  the 

_  _  fi 

curve  r£(x)  is  plotted  for  Mach  numbers  equal  to  8.5,  16,  and  22. 

The  pressure  has  been  kept  constant,  and  the  constant  of  Integration 

of  Eq.  (32)  has  been  set  equal  to  zero.  Mote  that  the  Influence  of 

the  cooling  by  expansion  Is  stronger  for  the  high  Mach  nuiri>ers  (and 

is  present  also  over  longer  distances  downstream). 

In  Appendix  E,  Figs.  E-1  through  E-6  give  the  function  r£(x)  for 

several  Mach  numbers  and  values  of  the  paraiseter  R,,  /C  . 

11  D££ 
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Fig.  5  —  The  growth  of  the  turbulent  core  assuming  constont  pressure  for  Cq  =  I.O,  M-16, 

and  several  R^,./C.  values 
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III.  AMALYSIS  TOR  THE  CASE  OF  CONES  AND  SLENDER  BODIES 

The  basic  difference  between  the  case  of  the  turbulent  wake  be* 
hind  a  blunt  body,  such  as  analyzed  In  the  previous  section,  and  that 
of  a  slender  one  is  that  In  the  first  geometry,  high  temperatures  are 
developed  at  the  stagnation  point  where  the  flow  comes  to  rest  by 
hitting  the  wall  perpendicularly,  whereas  In  the  second,  high  temper¬ 
atures  are  developed  In  the  boundary  layer  where  the  flow  comes  to 
rest  through  the  action  of  the  viscous  forces  which  act  parallel  to 
the  wall.  In  the  case  of  blunt  bodies  It  was  found  that  most  of  the 
momentum  deficiency  was  located  In  the  Invlscld  flow  through  the  ac¬ 
tion  of  the  shock  wave,  with  high  characteristic  enthalpies  prevail¬ 
ing  at  the  cross  section  x  ■■  equal  to  about  1/3  the  stagnation  en¬ 
thalpy  (H  (Y  -  1)1^ /6  ~  H^^/3).  The  distribution  of  the  drag  at 
the  neck  was  such  that  .  Consider  now  as  an  extreme  the 

case  of  a  thin  cone.  Even  at  high  Mach  numbers  the  oblique  shock 
wave  formed  at  the  tip  decelerates  the  flow  very  little,  and  hence 
the  Invlscld  enthalpies  are  very  low  compared  with  the  enthalpies 
prevailing  Inside  the  boundary  layer  or  at  the  neck. 

Figure  7  shows  schematically  the  relative  position  at  the  neck 
of  the  enthalpy  profiles  for  the  viscid  and  Invlscld  parts  of  the  two 
cases  under  discussion.  The  Invlscld  solutions  are  for  an  Infinite 
thin  cone  and  a  hemisphere -cylinder  and  are  modified  at  the  neck  by 
two  parabolic  profiles.  The  case  of  a  blunted  cone  will  lie  some¬ 
where  between  the  curves  for  a  pure  cone  and  hemisphere -cylinder.  As 
pointed  out  In  Ref.  9,  the  shock  %rave  in  the  neighborhood  of  the  tip 
Is  given  reasonably  well  by  the  cylindrical  blast-wave  theory  based 
on  the  curvature  of  the  nose,  whereas  far  from  the  body  the  same  is 
true  If  the  radius  of  the  base  is  taken  Into  consideration.  The  re¬ 
gion  In  between  Is  approximated  by  a  conical  shock  %dtich  allows  for 
a  slight  recompresslblllty  zone.  The  exact  shape  of  the  shock  wave 
Is  needed  for  the  determination  of  the  enthalpy  profile  at  the  cross 
section  X  •  where  the  pressure  Is  almost  ambient.  In  the  fashion 
of  Ref.  1  this  distribution  can  be  obtained  In  terms  of  the  shock 
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hoo 

.  7  —  Schematic  enthalpy  distribution  at  the  neck  in 
the  wake  of  a  blunt  body  and  a  thin  cone 
for  the  same  Mach  number 
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angle  3  by  using  the  regular  condltlens  across  the  shock  wave  and  ex¬ 
panding  isentropically  along  each  different -entropy  streamline.  This 
procedure  yields 

h  (R) 

Jo _ 

where 

sin^P  -  [l  +  (dR^/d^)"^]  (40) 

The  term  is  an  effective  specific-heat  ratio  which  accounts  for 

real  gas  effects.  For  blunt  bodies,  Y^  was  found  to  be  approximately 

equal  to  1.2  in  Ref.  1.  For  thin,  slender  cones,  Yj^  should  approach 

the  value  1.4,  since  the  temperatures  for  this  case  are  lower  by  far. 

The  distance  x  ■  x  is  located  where  the  angle  of  the  shock  wave  be- 
0 

comes  the  Mach  angle,  or 

*0  ■  rT±r=  ■  ^  ““ 

(dR^/dx/ 


Equation  (41)  is  an  exact  one,  valid  for  both  thin  and  blunt  bodies 

* 

but  not  for  infinite  cones. 

We  now  write  the  energy  equation  (Eq.  (13))  without  making  spe¬ 
cific  substitutions  for  the  functions  h^  and  A.  We  have 


dR. 


hf  +  A)  ^ 


dR,  - 


-8<dx 


f(^^7-^)(Y  -  1)M^ 


(42) 


The  term  h^  can  be  obtained  from  Eq.  (39)  by  making  R  ■  R^.  On 

A  finite  cone  will  have  a  shock  wave  bent  towards  the  center 
line  by  virtue  of  the  interaction  with  the  expansion  fan. 
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th«  other  hand,  the  expression  for  A(R^)  found  In  Appendix  C  needs 
no  modlflceClen  except  that  now  the  approximation  Is  net 

valid,  end  as  e  consequence  Eq.  (C-6)  must  be  used.  This  equation 
Is  reproduced  below  for  continuity  after  combination  with  Eq.  (29): 


A 


Cp^(Y  -  1)M^ 


(V  - 


(43) 


Equation  (42)  can  now  be  Integrated  by  using  Eqs.  (39)  and  (40)  for 

the  determination  of  h,,  and  Eq.  (43)  for  the  function  A. 

*  2  2 
Let  us  now  consider  as  an  example  the  case  for  which  M  sin  6  1* 

a  large  ntanber  In  order  for  the  following  approximation  of  Eq.  (39) 

to  be  permissible: (see  Appendix  F  for  a  detailed  discussion  of  this 

assumption) 


We  may  neglect  for  convenience  the  quantity  In  the  second  parenthesis, 
so  that 


h  (R) 

X 


1 


1  +  (dR^/dx)’^J 


(44) 


As  In  the  case  of  the  blunt  body  let  us  set 
h  (R) 

X 


(45) 
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The  prime  indicates  differentiation  with  respect  to  the  co-ordinate 
of  the  shock  wave  x,  not  to  be  confused  with  the  distance  downstream 
from  the  neck  over  which  we  study  the  growth  of  the  turbulent  core. 

Differentiation  of  Eqs.  (43)  and  (45)  and  substitution  in  Eq.  (42) 
yield,  after  some  algebra 


From  the  above  we  see  clearly  the  strong  dependence  of  the  growth  of 
the  turbulent  core  on  the  second  derivative  (or  curvature)  of  the 
shape  of  the  shock  wave. 

As  a  limiting  case  we  further  study  in  some  detail  the  case  of  a 
pure  cone  in  the  approximation  M^sin^^  »  1.  In  this  case  we  set  the 
factor  f  1,  and  from  Eq.  (4)  we  see  that  is  exactly  equal  to  the 

Howarth  variable  R.  Also,  If  0  is  the  angle  of  the  cone  (assumed 
small),  we  have  for  the  shape  of  the  shock 

As  •  »  "•  (  *  (^7) 

From  Eq.  (44)  we  also  find 
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H 


2/Yl 


(48) 


Equation  (42)  becomes 


The  second  expression  under  the  radical  on  the  left-hand  side,  repre¬ 
senting  the  inviscid  effects,  is  negligibly  small  for  small  0  compared 
with  the  third  part,  which  contains  the  strong  viscous  effects.  We 
now  Integrate  Eq.  (49)  in  the  limiting  case  of  small  x.  We  distin¬ 
guish  the  following  cases: 

A.  Let  us  assume  that  Cjj  * 

This  is  the  case  of  very  thin  cones  (a  few  degrees  half 
angle)  and  small  Reynolds  numbers.  Pressure  variation 
apart,  Eq.  (49)  yields  after  integration 


(50) 


B.  Let  C  «  C_  «  l.O. 
fi  ** 

In  this  case  the  denominator  of  the  last  part  under  the 

radical  is  approximated  by  R,.^/C  .  After  integration 

we  have  ^i 
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^(Y  -  DM^C 


(x  -  x^)  + 


(51) 


For  completeness  and  for  the  benefit  of  ready  comparison  we  give  be¬ 
low  the  equivalent  result  for  blunt  bodies  (Eq.  (36)).  We  assxune 
(1  +  f)/2f  ^1.0,  but  ve  Incorporate  the  pressure  term  In  an  average 
sense,  since  for  blunt  bodies,  especially  at  high  Mach  numbers.  It  Is 
significant. 


-  -  hil 
% 


2/3 


(52) 


The  decay  of  the  enthalpy  along  the  axis  Is  given  by  Eq.  (19)  as 
follows: 


hM 


2/Y, 


■Cy  - 1)^ 


(53) 


The  equivalent  result  for  a  blunt  body  Is 


liM 


(Y  -  1)M^ 


(54) 


From  Eq.  (S3)  we  see  that  in  the  case  of  a  pure  cone  the  part  of  the 
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enthalpy  due  to  the  shock  wave  is  negligibly  small  compared  to  the 

part  due  to  the  viscous  effecta.  From  Bq>  (54)  we  see  that  for  small 

X  the  reverse  Is  true  for  blunt  bodies. 

2 

We  now  examine  the  quantity  R,.  /C-  In  the  case  of  cones.  From 

ri 

Bq.  (34),  valid  for  all  cases,  we  note  that  h^(e)  »  h,,,  and  hence 
Eq.  (34)  reduces  In  terms  of  the  parameter  \  to 


""^fi 


i 

\ 


(55) 


Bearing  the  above  In  mind,  we  use  Eqs.  (50),  (51),  and  (55)  In  Eq.  (53), 
and  Eqs.  (35)  and  (52)  In  Eq.  (54);  we  find  the  following: 

For  the  case  of  a  cone  with  -  C.  «  1  and  for  small  x 

D  » .  U 


h(o) 

h. 


■(1.- 


8/C 


(x  -  x^) 


(56) 


For  the  case  of  a  cone  with  «Cj,«  1  and  for  small  x 


For  the  case  of  a  blunt  body  with  Cp  « 


fi 


1.0 
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(Y  -  uid. 
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72(3X  -  ir/< 


(Y  -  DM^Cjj  L 


2I0I 


javg 


2/3 


(58) 


The  following  conclusions  can  now  be  drawn  from  the  above  equations: 

1.  The  level  of  the  initial  enthalpy  at  the  neck  Is  of  the  same 
order  for  both  blunt  and  slender  bodies,  since  at  the  neck  there  is  a 
tendency  for  the  stagnation  enthalpy  to  be  recovered  at  the  point  where 
the  dividing  stream  line  of  the  free  shear  layer  crosses  the  axis. 

2.  In  the  cases  of  cones  or  blunt  bodies,  the  smaller  the  amount 
of  total  drag,  the  higher  the  decay  of  the  enthalpy,  as  one  would  ex¬ 
pect  Intuitively.  This  Is  shown  by  the  fact  that  the  factor  which 
multiplies  (x  -  x^)  in  Eqs.  (56),  (57),  and  (58)  increases  with  smal¬ 
ler  Cj,. 

3.  Comparison  of  Eqs.  (57)  and  (58),  representing  the  most  usual 

cases  In  practice,  shows  that  for  the  same  M,  and  C_  the  enthalpy 

fi 

decays  faster  for  a  cone  with  «  1.0.  Examination  of  the  quanti¬ 
ties  multiplying  (x  -  x^)  in  Eqs.  (57)  and  (58)  shows  that  their  ratio 
(for  a  value  \  •  1/2)  Is  equal  to 


where  the  subscripts  b  and  c  stand  for  the  blunt  and  conical  bodies; 
the  pressure  term  refers  to  the  blunt  case,  since  for  thin  cones,  cool¬ 
ing  by  expansion  is  not  significant.  One  can  see  that  the  pressure 
cooling  for  the  blunt  bodies  competes  with  the  reduced  drag  mechanism 
for  slender  bodies,  especially  at  the  high  Hach  ntimbers.  From  Ref.  4 
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at  M  ■  22  we  estimate  ““ 

(Cjj)^  ■  0.1  the  quantity  multiplying  (x  -  in  the  conical  case  will 

be  a  little  higher  when  compared  to  the  blunt  case.  On  the  other  hand 
one  can  see  again  by  observing  Eqs.  (57)  and  (58)  that  the  expression 
containing  x  In  the  braces  In  Eq.  (57)  dominates  faster  the  number  to 
Its  left,  since  this  number  Is  smaller  than  6  (which  corresponds  to  the 
blunt  case)  as  long  as  X  >  1/3.  The  above  effects  are  shown  qualita¬ 
tively  in  Fig.  8. 

For  the  case  of  a  blunted  cone,  the  Invlscld  enthalpy  distribution 

will  look  qualitatively  as  shown  In  Fig.  9,  due  to  the  fact  that  the 

(9) 

flow  suffers  a  recorapresslon  In  the  conical  part  of  the  shock.  Since 
Inflection  points  Influence  the  growth  of  the  turbulent  core  (see  Eq. 
(46)),  through  the  second  derivative  R^',  the  growth  curve  r^(x)  will 
contain  a  number  of  wiggles  depending  on  the  relative  position  of  the 
"viscous"  parabolic  profile  with  respect  to  the  Invlscld  one.  It  might 
be  Instructive  to  compare  a  spherically  blunted  cone  and  a  pure  cone 
having  the  same  total  drag  coefficient  and  flying  at  the  same  altitude 
and  Mach  number.  For  the  reasons  explained  above,  we  expect  the  growth 
curve  r^(x)  of  the  blunted  cone  to  exhibit  Irregularities  due  to  over- 
expansions  and  recompressions.  Eventually,  at  high  distances  x,  when 
the  wake  behaves  like  an  incompressible  one,  both  curves  will  tend  asymp¬ 
totically  toward  the  same  straight  line  In  the  log-log  representation. 
This  effect  Is  shown  In  Fig.  10.  On  the  other  hand,  It  Is  easy  to  see 
from  Eq.  (52)  that  the  decay  of  the  enthalpy  at  the  axis  will  be  slower 
for  the  blunted  cone  In  the  Initial  stages  than  for  the  pure  cone,  re¬ 
suming  the  same  rate  later  on  when  the  compressibility  effects  are  again 
absent.  Numerical  solutions  for  several  conical  geometries  are  reported 
In  Ref.  12 . 


This  has  been  tacitly  assumed  In  writing  Eq.  (35);  physically  It 
means  that  the  enthalpy  at  the  neck  must  be  higher  than  the  enthalpy  at 
the  cross  section  where  the  pressure  Is  ambient. 


X 


Fig.  8  —  Schematic  decay  of  the  enthalpy  at  the  axis 
of  a  wake  behind  a  blunt  body  and  a  cone  for 
the  same  Mach  number 


Fig.  9  —  Schematic  enthalpy  profile  at  the  neck 
In  the  wake  of  a  blunted  cone 
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Blunted  cone 


FIg.lO — Schematic  shape  and  relative  growth  of  the 
turbulent  core  behind  a  blunted  cone  and  a 
cone  of  the  same  total  drag  coefficient 
for  the  same  Mach  number 
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IV.  CONCLUSIOHS 

An  attempt  was  made  to  study  the  growth  of  the  turbulent  core  in 
the  wake  behind  blunt  and  slender  bodies  flying  at  hypersonic  speeds. 
An  exact  analysis  Is  not  possible  because  of  Incomplete  understanding 
of  the  following:  (a)  the  flow  In  the  near  wake  and  the  neck  where 
a  recompression  occurs  through  the  trailing  shock;  (b)  the  mechanism 
of  transition  to  turbulence;  (c)  the  structure  of  fully  established 
turbulent  flow  In  a  compressible  wake. 

In  order  to  compensate  for  problem  (a)  above,  we  treat  the  quan¬ 
tity  of  the  enthalpy  at  the  center  of  the  neck  as  an  additional  param¬ 
eter,  and  we  proceed  to  study  Its  effect.  Problem  (b)  Is  more  serious. 

In  the  present  work  we  assume  that  fully  established  turbulence  starts 

(10) 

at  the  neck.  Evidence  from  ballistic-range  experiments  suggests 
that  this  Is  the  case  at  the  higher  densities,  whereas  at  lower  den¬ 
sities  the  transition  point  (roughly  defined  by  the  distance  where  the 
Intermlttency  region  appears)  could  appear  many  diameters  downstream. 
The  theory  presented  here  can  easily  be  modified  to  take  Into  account 
this  last  effect,  provided  that  the  Initial  conditions  from  the  point 
where  the  turbulent  flow  Is  fully  established  are  known.  For  problem 
(c) ,  the  best  that  can  be  done  Is  to  extend  our  semlemplrlcal  knowledge 
on  the  subsonic  wakes  to  the  hypersonic  case  as  suggested  In  Ref.  3. 
Only  detailed  comparison  of  many  of  the  quantities  predicted  by  the 
theory  will  show  to  what  extent  this  procedure  Is  correct.  It  must 
be  admitted  that  the  growth  of  the  turbulent  core  Is  not  the  most  sen¬ 
sitive  Indicator  (temperature  decay  would  be  better) ,  but  It  appears 
that  the  theory  Is  In  good  agreement  with  ballistic-range  experiments. 

In  this  work,  a  number  of  approximations  have  been  made  so  that 
the  final  results  would  be  In  a  closed  form.  In  order  that  the  differ¬ 
ent  trends  for  a  great  variety  of  cases  could  be  appreciated  without 
the  need  for  detailed  numerical  calculations.  On  the  basis  of  such  a 
study,  the  following  conclusions  can  be  drawn: 

1.  All  bodies  produce  rates  of  growth  for  the  turbulent  core 
which  are  higher  for  higher  Hach  numbers  and  hotter  viscous  cores. 

The  effects  of  both  Mach  number  and  enthalpies  prevailing  at  the  neck 
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•re  cancelled  out  •  few  hundred  diameters  downstream,  where  the  core 
grows  with  the  1/3  power  of  distance  multiplied  by  the  drag  coefficient 
and  the  cross  section. 

2.  Blunted  cones  produce  growths  with  more  Irregularities  be¬ 
cause  this  geometry  creates  overexpansions  and  recompressions. 

3.  The  level  of  enthalpies  at  the  neck  for  all  bodies  flying  at 
the  same  altitude  and  Mach  number  Is  of  the  same  order  of  magnitude 
(smaller  for  slender  bodies  after  allowance  Is  made  for  boundary-layer 
effects) ,  but  the  decay  of  the  enthalpy  along  the  axis  of  the  wake  Is 
faster  for  a  cone  than  for  a  blunted  body,  even  when  both  cone  and 
blunted  body  may  have  the  same  total  drag  coefficient. 
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Appendix  A 

CALCULATIOM  OP  THE  VELOCITY  PREVAILING  AT  THE  POINT  ON  THE 
AXIS  WHERE  THE  PRESSURE  IS  AUCST  AMBIENT  IN 
THE  IMVISCID  MAKE  OP  A  BLUMT  BODY 

In  Ref.  1  It  was  shovm  that  the  value  of  the  enthalpy  on  the 

axis  of  the  wake  of  a  blunt  body  and  at  the  distance  x  where  the 

o 

pressure  is  ambient  is  given  as 

‘'x  /  \  p 

Jo _  ^ 

h,  \  2  /  1_  2Y  _ 

It  was  also  shown  that(  independent  of  Mach  number,  1.2.  It  is 

easy  to  see  that  the  quantity  in  brackets  in  Bq.  (A-1)  is  very  close 
to  one .  Hence 


1-  fr 


'L/  .  1/Y, 


(M^) 


(A-1) 


h,  (0) 

Ji _ 


(A-2) 


Assuming  conservation  of  the  stagnation  enthalpy,  we 


can  show  that 


(1) 


f 


2 


1 


2 

(Y  -  1)»^ 


h  (0) 

*0 


1 


(A-3) 


Comparing  the  last  t%ro  equations,  we  find 


f 


(A-4) 
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This  relation  shows  that  £  is  highly  Insensitive  to  variations  In 
Mach  number.  Some  values  are  computed  In  the  following  table: 


M  i 

10 . 0.73 

15 . 0.77 

20 . 0.80 

25 . 0.81 

30 . 0.82 

35 . 0.83 
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Appendix  B 

INFLUENCE  OF  THE  LAMINAR  COOLING  AT  THE  FRONT 
OF  THE  TURBULENT  WAKE 


From  Ref.  1  it  Is  found  that 


hfo.Xl 

h(o,o) 


(■ 


(B-1) 


where 


X 


Re 


with  Re 


P»Vo 


(B-2) 


Let  us  choose 

~  2, 


f  =■  0.8,  Pr  -  0.7,  =*•  Cp 


Then  the  above  becomes 


h(o,o) 


(B-3) 


_x2  2  -1 

Using  the  approximation  e  (1  x  )  ,  which  Is  compatible  with 

Eq.  (10)  and  the  findings  of  Ref.  1, 


The  last  term  In  the  denominator  can  be  neglected  If 


ii  2L. 
Re 


«  R‘ 


(B-5) 


Now  at  large  x  we  have  x  2R 


3 


(see  Fig.  3);  hence  the  above  yields 


or 


2i2^ 

Re 


(B-6) 


From  the  data  of  Ref*  8  the  minimum  Reynolds  number  of  all  ex* 

A 

perlments  reported  Is  approximately  equal  to  10  .  Setting  ~ 
we  find 


X  « 


(B-7) 


This  Is  Indeed  a  very  high  number 
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Appendix  C 

EVALUATION  OF  THE  FUNCTIONS  A  AND  C-  IN 
TERMS  OF  THE  HOWARTH  VARIABLE  R 


The  total  drag  at  any  given  distance  downstream  X|  associated 

♦ 

with  the  turbulent  core  only,  is 


(Drag)  ■  2n  J  ^  pu(u  -  u  )  rdr 


(C-1) 


Through  Eqs.  (16),  (19),  and  (20),  the  above  becomes 


(Drag) 


^.A(R)r^(R^^)^  Pff 


m 


(C-2) 


We  define  a  local  drag  coefficient  as 


(Drag) 


TTr 


(C-3) 


Equation  (C-2),  after  solving  for  A(R)  and  using  Eq.  (C-3),  becomes 

.2 


A(R)  - 


Cp  PJJ. 


Since  A(R)  will  only  be  used  in  the  neighborhood  of  the  axis,  we  have 
from  Eq.  (25) 


* 

Assuming  uniform  pressure  everywhere  and  zero  velocity  com¬ 
ponent  in  the  radial  direction. 


-2  \  2f  ; 


(C-4) 


For  simplicity  we  assume  from  now  on  that  (1  +  f)/2f  “■  1.0. 

Let  us  assume,  for  the  purposes  of  this  Appendix  only,  that  the 
enthalpy  distribution,  which  includes  an  equivalent  drag  coefficient 
(C_  )  within  the  distance  R,, ,  is  given  by  the  relation 

Dfi 


We  set 


h(R)  .  -RVe 
h(R  -  o) 


(C-5) 


It  Is  obvious  that  as  R^  **  *»  "*  Cp  as  it  should.  Integration 

yields  ^ 


-  1  -  e  * 


For  small  R^  we  have 


1  +  ^  *  »f 


It  2 

Mote  that  the  approximation  e  "  1/(1  +  x  )  Is^good  for  almost 

-x^  2 

all  values  of  x  between  zero  and  infinity,  vihereas  e  1  •  x  ob¬ 
viously  Is  limited  to  values  of  x  that  are  very  small  compared  to  one. 
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For  blunt  bodies  for  which  ,  the  above  simplifies  to 


(C-7) 
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Appendix  D 

CALCULATIOM  OF  r  IN  TERMS  OF  THE  HDMARTH  VARIABLE  R 
TOR  BLUMT  BODIES 


The  Howarth  variable  R  is  defined  as 


p  RdR  ■  Prdr 


Assuming  an  ideal  gas,  the  above  is  written  as 


“  h  .-2  -2 

—  7”  dR  "dr 
P  K 


Through  Eq.  (10),  this  reduces  to 


.-2 


H 


dr  “  ■  —  /I  +  , 

P  I 

1  +  ~ 


dR- 


If  we  assume  the  pressure  p  to  be  constant  for  all  R,  the  above  in¬ 
tegrates  to 
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Appendix  E 

STUDY  OF  THE  FUNCTION  r^(x)  FOR  DIFFERENT  M4CH  NUMBERS  AND 

DIFFERENT  ENTHALPIES  FREVAttlNG  AT  THE 

CENTER  OF  THE  NECK  OF  A  BLUNT  BODY 

Figures  E-1  through  E-6  show  the  function  r^Cx)  for  Mach  numbers 

5y  7,  10,  13,  19,  and  25.  For  each  Mach  number,  the  value  of  the 
2 

parameter  K,JC-  ,  which  depends  only  on  the  percentage  \  of  the  stag- 
ri 

nation  enthalpy  prevailing  at  the  center  of  the  neck  according  to  Eq. 
(35),  has  been  set  equal  to  3,  4,  6,  8,  10,  and  12. 
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Fig.E-2 — The  growth  of  the  turbulent  core  assuming  constant  pressure 
for  Cq=I.O,  M=7,  and  several  R^fj/c^  values 


50 


Fig.E-3 — ^The  growth  of  the  turbulent  core  assuming  constant  pressure 
for  Cq=I.O,  M>IO,and  several  R^fj/Cp^.  values 
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Fig.E-4 — The  growth  of  the  turbulent  core  assuming  constant  pressure 
for  Cq=I.O,  M=l3,and  several  R^fj/Cp  values 
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Appendix  F 

APPROXIMATION  OF  EQ.  (39)  WITH  EO.  f44) 


We  need  only  to  make 


M^sin^e  » 

because  we  satisfy  simultaneously  the  condition 

M^sln^P  » 

The  above  statement  follows  from  the  Inequality 


(F-1) 


(F-2) 


(F-3) 


which  holds  If 


3  y  <  3  (F-4) 

The  above  Is  always  true. 

Romlg  In  Ref.  11  has  found  for  conical  flows  allowing  for  dis¬ 
sociation  the  approximation 

M  8ln3  “•  0.39  +  1.03  M  sin  Qf  (F-5) 


In  the  above,  M  Is  a  reference  Mach  number  defined  by  the  relation 


M 


(F-6) 


with  Tj^  ■  491. 7°K,  and  ®  the  angle  between  the  cone  axis  and  the  cone 
surface. 
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Making  ^  and  using  Eq.  (F-5)  In  (F-1)  we  find  after  some 
algebra  that  Inequality  (P-l)  Is  equivalent  to 

M  sin  Of  »  7-^*^  -  0.38  (P-7) 

V  Y  -  1 

for  Y  "  1>A  the  above  becomes 

M  sin  Of  »  1.80  (F-8) 

As  an  example,  for  a  20-deg  cone  M  »  5.25,  and  for  a  10-deg  cone 
M  »  10.4. 
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